The influence of magnetic flelds on the thermodynamics of 

primordial star formation 



Dominik R. G. Schleicher^, Daniele GaUi^, Simon C. O. Glover^, Robi Banerjee^, 
Francesco Paha ^, Raffaeha Schneider^, Rah S. Klessen^ 

dschleic@ita.uni-heidelberg.de , galliOarcetri . astro . it , sglover@ita.uni-heidelberg.de , 
baiierjee@ita.uni-heidelberg.de , palla@arcetri . astro . it , raf f afiarcetri . astro . it , 

rklessenSita . uni-heidelberg . de 



ABSTRACT 



We explore the effects of magnetic energy dissipation on the formation of the first stars. For 
this purpose, we follow the evolution of primordial chemistry in the presence of magnetic fields in 
the post-recombination universe until the formation of the first virialized halos. From the point 
of virialization, we follow the protostellar collapse up to densities of ~ 10^^ cm~^ in a one-zone 
model. In the intergalactic medium (IGM), comoving field strengths of > 0.1 nG lead to Jeans 
masses of 10^ Mq or more and thus delay gravitational collapse in the first halos until they 
are sufficiently massive. During protostellar collapse, we find that the temperature minimum at 
densities of ~ 10^ cm~'^ does not change significantly, such that the characteristic mass scale 
for fragmentation is not affected. However, we find a significant temperature increase at higher 
densities for comoving field strengths of > 0.1 nG. This may delay gravitational collapse, in 
particular at densities of 10^ cm~^, where the proton abundance drops rapidly and the main 
contribution to the ambipolar diffusion resistivity is due to collisions with Li+. We further 
explore how the thermal evolution depends on the scaling relation of magnetic field strength with 
density. While the effects are already significant for our fiducial model with B esc p0-5-0-57 ^ ^^le 
temperature may increase even further for steeper relations and lead to the complete dissociation 
of H2 at densities of ~ 10^^ cm~'^ for a scaling with B oc p^'^. The correct modeling of this 
relation is therefore very important, as the increase in temperature enhances the subsequent 
accretion rate onto the protostar. Our model confirms that initial weak magnetic fields may 
be amplified considerably during gravitational collapse and become dynamically relevant. For 
instance, a comoving field strength above 10~^ nG will be amplified above the critical value for 
the onset of jets which can magnetize the IGM. 
Introduction 



In the absence of magnetic fields, protostellar 
collapse during first star formation is well under- 
stood and ha s been followed in detail by numerical 
simul ations (jAbel et al.l 120021 : iBromm fc Larson 
20041 ). The main physical processes, includ- 



tation and protoste llar feedback hav e been dis- 
cussed in detail by ICiardi fc Ferraral (jiooi) and 
Gloveii (j2005l ). Recent numerical simulations 



ing chemistry, cooling, the extent of fragmen- 
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confirme d this scenario with unprecedented ac- 
curacy IjYoshida et al.l 120081 ). The first stars 
are expected to be very massive (with charac- 
teristic masses of ~ 100 Mq), consistent with 
constraints from the observed reionization opti- 
cal de pth (jWvithe fc Loebll2003t ISchleicher etal 



2008aD . 



However, it may be important to consider the 
effects of put ative magneti c fields. Recent ob- 
servations by iBernet et al. showed that 
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strong magnetic fields {B ^ were present 

in normal galaxies at z '--^ 3, where relatively lit- 
tle time was available for a dynamo to operate. 
In our galaxy, the field strength is 3 — 4 fiG, 
and it is coherent over kpc scales, with alternat- 
ing directions in the arm and interarm regions 
(e.g. iKronbergl Il994l : iHanl l2008l ). Such alterna- 
tions are expected for magnetic fields of primordial 
origin, but m ore difficult to explain in m ost dy- 
namo models ( Grasso fc Rubinstein 200 More- 
over, it is not clear whether large-scale dynamos 
are efficient, as the small-scale magnetic fields 
are produced at a faster rate and lead to satu- 
rati on before a significa n t large-scale fi eld builds 
up (|Kulsrud et al.lll997l ). IZweibell (|2006l ) confirms 
problems with explaining the observed magnetic 
field strength from dynamo theory in our galaxy, 
while models based on stellar magnetic fields have 
difficulties with explaining the large-scale coher- 
ence. On the other hand, magnetic fields ob- 
served in some other spiral galaxies appear to be in 
agreement w ith the predictions of dynamo theory 
(lBeckll2009h . The main arguments against pri- 
mordial magnetic f i elds h ave been considered by 
Kulsrud k Zweibell (|2008h . who find that they are 
too uncertain to rule out this possibility. Strong 
magnetic fields have also been detected in the po- 
larized disk in M 31 ( Berkhuijsen et aT]l2003h . 

At redshifts higher than z > 3, it is only 
possible to derive upper limits on the magnetic 
field strength. Observations of small-scale cos- 
mic microwave background (CMB) anisotropy 
yield an upper li mit of 4.7 nG to the comov- 
ing field strength ( Yamazaki et al.l l2006l ). Ad- 
ditional constraints can be derived from the mea- 
surement of tJgj which describes the root-mean- 
square of the density fluctuatio ns at 8h~^ comov- 



1978HKim et al.lll996l :| Sethi fc Subramanianllm 



ing Mpc (jYamazaki et al.ll2008l ). The inferred field 
strength is similar, but depends on the assumed 
shape of the power spectrum. Primordial nucle- 
osynthesis constrains the co moving field strength 
to b e less than ~ 1 (jGrasso fc Rubinstein 
Il996l) . while reionization yields upper limits of 
0.7— 3 nG, depending on assumptions on the stel- 
lar population that is responsible for rcionizing 
the universe ISchleicher et al.1 l2008al ) . As shown 
in previous works, these upper limits are rela- 
tively weak, but may imply a significant impact 
on the thermal evolution in the pos trecombina- 
tion universe and structure formation ( WassermanI 



iTashiro fc Sueivama' '20 061: iTashiro et al""" '2006 
^Scthi ct al. 2008; Sc hleicher et all l200 8a. 200S 
Numerical simulations based on an SPH-MHD 
version of Gadget indicate that a comoving field 
strength of 2 x 10"'^ nG might be sufficient to 
reprodu ce the observed magnet ic fields in galaxy 
clusters ( Dolag et al.ll2002 . 20051 ). However, larger 
initial field strengths also appear to be consistent 
with the observed data, whic h shows a large de - 
gree of scatter (see Fig. 10 in lDolag et al.ll2002l) . 
Because these results may also change when in- 
cluding processes like AGN feedback, cooling and 
heat conduction, more work is required before be- 
ing able to draw definite conclusions about the 
upper limit to the field strength. 

For the early universe, a number of viable sce- 
narios exist that can produce primordial fields. 
Mechanisms to generate such fi elds in the epoch of 
i nflati on have been proposed bv lTurner fc Widrow 
(jl988| ). and more recent works confirm that strong 
magnetic fields of up to 1 n G (comoving) may 
be p roduce d in this epoch (Bertolami fc Mot 
1999.: Bamba et al.l 120081 : [Cam panelli ct al. 20a 



Campanellill2008f) . These models require that con- 



formal invariance is broken explicitly or implicitly 
(jParker 1968 ). Magnetic fields may al so form dur- 



ing t he electroweak phase transition (jBavm et al 
19961 ). especially if it is a phase-transition of first 



order, as requi red in the electroweak baryoge- 
nesis scenario (|Riotto fc TroddenI Il999l ). simi- 
larly, the QCD phase transiti on may give rise to 



strong primordial fiel ds (Ou ashnock et al. 119891 : 



[Ch eng fc 01intol[l99l . Particularly strong fields 
(~ 1 nG, comoving) can be generated if hydrody- 
namical instabilities were present at this epoch, 
as ex pected for a lar ge range of model parame- 



ters (|Sigl et al.l 119971 ). The further evolution of 



magnetic fields generated at the QCD or elec- 
troweak phase transit i on ha s been calculated by 



Banerjee fc Jedamzild (|2004[ ). who find that the 



coherence length may increase further if the mag- 
netic helicity is non-zero. 

There is thus a reasonable possibility that mag- 
netic fields were present in the universe after re- 
combination and during the formation of the first 
stars. Previous works have started to explore the 
effects of magnetic fields, both those of primordial 
origin, as well as those created du ring structure 
formation or in the protostellar disk. iMaki fc Susa 
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(|2004l 120071) investigated the coupling of the mag- 
netic field to the primordial gas during collapse, 
finding that the field is frozen to the gas and may 
drive the magnetor otational instability (MR I) in 
the accretion disk. Tan k, BlackmanI (|2004l ) pro- 
posed that gravitational instabilities in the accre- 
tion disk may give rise to turbulence and drive 
a dynamo that increases t he field streng t h un- 
til equipartition is reached. Silk fc Langer ( 2006h 
discussed whether the MRI can drive a dynamo 
in the protostellar disk, suggesting that magnetic 
feedback may give rise to a distribution of stellar 
masses closer to the present-day initial mass func- 
tion (IMF) . The effects of such fee dback have been 
explored by iMachida et al. I (|200d) . who find that 
jets may blow away up to 10% of the accr eting 
matter. On the other hand, IXu et al. I (|2008l) cal- 
culated the generation of magn etic fields by the 
Biermann battery mechanism ( Biermamil 1950l ) 
during the formation of the first stars, and find 
that such fields are not strong enough to affect 
the dynamics of first star formation. Simula- 
tions of present-day star formation indicate that 
magnetic fields suppress fr agmentation of cluster- 
form in g molecular c ores (|Hennebelle fc Tevssier 



20081: iPrice fc Batd 120081 ) and extract a con 



siderable amount of angula r momentum from 
the protostellar disk (e.g., iBanerjee fc Pudritz 



20071 ) helping to form more massive stars. Mas- 
sive stars may play an important role in creat- 
ing and spreading magn etic fields on timescales 
of a f ew million years dBisnovatvi-Kogan et al 



1973al lbl: iBisnovatvi-Kogan fc Lamzin 



recently shown by Hanasz et al.l 



19771) 



2009 V. 



As 
the 



cosmi c ray-driven dynamo proposed by I Parker 
()l992l) could work on timescales faster than 
the rotation timescale of a galaxy, and there- 
fore be important even for the earliest gala xies 
(|Bouwens fc mingworthll2006l : live et al.ll200d ). 

The effects of magnetic fields on the first stars 
are therefore still subject to significant uncertain- 
ties. In this work, we point out that magnetic 
fields may not only directly affect the dynamics 
by launching jets or driving the MRI in the accre- 
tion disk, as previously suggested, but also indi- 
rectly through changes in the chemistry of the gas 
and additional heat input from magnetic energy 
dissipation. We focus on comoving field strengths 
of up to ~ 1 nG, as stronger fields are difficult 
to generate in the early universe and current con- 



straints would not allow for their existence. We 
calculate the evolution of primordial chemistry in 
the presence of magnetic fields until the forma- 
tion of virialized halos in § [21 In § [31 we follow 
the thermal evolution of the gas during collapse 
in a one-zone model and show the impact of the 
magnetic field on the thermal and magnetic Jeans 
mass. The implications of our results are discussed 
in§[4l 

2. From the large-scale IGM to virialized 
halos 

Here, we explore the effects of magnetic fields 
on the chemical evolution in the post-recombination 
universe until the formation of the first virialized 
halos at z ~ 20. This calculation provides the ini- 
tial conditions used for our models of protostellar 
collapse that is described in § [3l 



2.1. Model description 

We calculate the thermal evolution of the IGM 
in the postrecombination universe with a modi- 
fied v ersion of the RECFAST codeQ (jSeager et al 



19991 ). a simplified but accurate version of a de- 
tailed code that follows hundreds of energy lev- 
els for H, He and He^ and self-consistently calcu- 



lates the background radiation field (jSeager et al 
200flf ). The r ecombinatio n calcu lation was recently 



updated bv IWong et al. ( 2008 ). and the physics 
of magnetic fields i n the IGM and of reion ization 
have been added bv lSchleicher et aL ( 2008al ) . Here 
we extend this model to follow the thermal evolu- 
tion until virialization and to self-consistently fol- 
low the chemical evolution until that point. At 
high redshift z > 40, the universe is close to ho- 
mogeneous, and the evolution of the temperature 
T is given as 



dz 



Xe (T - Trad) 



3i7(z)(l 
2T 



f z)meC 1 -I- /hc + Xe 
2(iAD — -^cool) 



1 + z ZnkBH{z){l + z)'' 



(1) 



where £ad is the heating term due to ambipo- 
lar diffusion (AD), icooi is the cooling term in- 
cluding Lyman a, bremsstrahlung and recom- 
bi nation cooling based on the cooling functions 
of lAnninos eTaD (Il997l) . In Eq. (|l]), ut is the 
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Thomson scattering cross section, the Stefan- 
Boltzmann radiation constant, me the electron 
mass, c the speed of hght, ks Boltzmann's con- 
stant, n the total number density, Xe = Ui^/nn 
the electron fraction per hydrogen atom. Trad the 
CMB temperature, H{z) is the Hubble factor and 
/hb is the number ratio of He and H nuclei. The 
latter can be obtained as /ho = ^/4(1 — Yp) from 
the mass fraction Yp of He with respect to the to- 
tal baryonic mass. T he AD heating rate is given 
as (|Pinto et al.ll2008f) 



^AD 



An 



(v X X B/B 



where 7jad is given as 



(2) 



(3) 



In this expression, the sum includes all neutral 
species tt-, and 77AD,n 

denotes the AD resistivity of 
the neutral species n. We calculate these resistivi- 
ties based on t he multifluid appro ach described in 
Appendix B of iPinto et al.l (|2008l ). This approach 
is a more genera l version of the formalism used by 
ISchleicher et al. ( 2008a ). and is particularly con- 
venient to calculate the AD heating rate in the 
high-density regime in § [3l 

Using this approach, we have calculated both 
the ambipolar diffusion heating rate as well as the 
Ohmic heating rate from the Hall parameters Pen , 
defined as 



f3sn = 



qsB 



iav 



(4) 



In this context, the index s denotes ionized species 
and n neutral species. Qs is therefore the charge 
of the ionized particle, rus its mass, to„ the mass 
of the neutral particle, pn the density of the neu- 
tral species, and ( 

^'^}sn the zero drift velocity mo- 
mentum transfer coefficients for these species (the 
zero drift approximation holds in the absence of 
shocks) . 

From these Hall parameters, it is possible to 
calculate the conductivity parallel to the electric 
field, (T||, as well as the Pedersen and Hall con- 
ductivities, dp and an, respectively. In general, 
these are calculated as sums over different pairs 
of ions and neutrals. In this paragraph, we focus 
only on the dominant contribution, which typi- 
cally comes from from collisions between the most 



abundant ion and neutral species. Then, the re- 
sistivities scale as cr|| oc /3s„, ap (x /3sn/(l + /3sn) 
and aH oc 1/(1 + /3s„). We further checked that 
(3sn >> 1 in the different regimes of our calcu- 
lation. Therefore, cr|| >> ap >> ch- For this 
reason, the ambipolar diffusion resistivity scales 
as rjAD (crp/(crp + a'jj) — l/c||) cx l/ap, while 
the Ohmic resistivity scales as rjo oc l/cr||. We 
therefore find that rjAD » Vo- We have checked 
this in more detail by evaluating the ambipolar 
diffusion and Ohmic heating rates during the cal- 
culation. As ambipolar diffusion always appeared 
to be the dominant magnetic energy dissipation 
mechanism in the cases considered here, we will 
focus on this contribution in the subsequent dis- 
cussion. Further details regarding th e mult ifluid 
approach can be found in lPinto et al.l (|2008r ). 

In the IGM, the dominant contributions to the 
total resistivity are the resistivities of atomic hy- 
drogen and helium due to collisions with protons. 
For their calculation, we adopt the mom entum 
transfer coefficients of Pinto &: Galli 

As the power spectrum of the magnetic field 
is unknown and Eq. ([2|) cannot be solved exactly, 
we adopt a simple and intuitive approach to esti- 
mate the differential operator for a given average 
magnetic field B with coherence length Lb- The 
heating rate can then be evaluated as 



VadB^ 



(5) 



In the IGM, one can show that ?7ad oc B^, so 
that w e recover the same depe ndence on B and 



Lb as ISchleicher et al.l (|2008a^ . The coherence 
length Lb is in principle a free parameter that 
depends on the generation mechanism of the mag- 
netic field and its subsequent evolution. It is 
constrained through the fact that tangled mag- 
netic fields are strongly damped by radiative vis- 
cosity in the pre-recombination universe on scales 
sm a.ller than the Alfven damping scale k~l^ given 



by ('jedamzik et al."l998': 'Subram anian Sz Barrow 
,1998 : .Scshadri fc Subr amanian 2o6l| 



234 Mpc"^ 



0.02 



2\ 1/2 



Bp 
1 nG 



0.7 



"(13 



1/4 



1/4 



(6) 



where Bq = B/{1 + z)^ denotes the comoving 
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magnetic field, ^Im and fli, are the cosmologi- 
cal density parameters for the total and baryonic 
mass, and h is the Hubble constant in units of 
100 km s~^ Mpc~^. In fact, we expect that fluc- 
tuations of the magnetic field may be present on 
larger scales as well. We thus estimate the heating 
rate by adopting Lb = k~l^ (l + z)~^, which yields 
the do minant contribution. Sethi fc SubramanianI 
also considered contributions from decay- 
ing MHD turbulence. However, these contribu- 
tions were found to be negligible compared to the 
AD heating terms, and highly uncertain, since 
only mild and sub-Alfvenic turbulence can be ex- 
pected at these early times. Even in the first star- 
forming halos, numerical hydrodynamics simula- 
tions show that the gas flow is laminar and only a 



small amount of turbulence is present (jAbel et al 
2OO2I : iBromm k Larsonll200l ). We therefore ne- 
glect this contribution here. 

The additional heat input provided by AD af- 
fects the evolution of the ionized fraction Xp of 
hydrogen, which is given as 



dxp 
dz 



H{z)(l + z)[l + KH{hH 
[1 + KH^Hnnil- Xp)] - 



- l3H)nH{l - Xp)] 



^riHX 



p 



H{z){i + zy 



(7) 



Here, tlh is the number density of hydrogen 
atoms and ions, 



hp Planck's constant, 



collisional ionization rate coefficient (Abel et al 



fcinn is th e 



1997). For the further details of notation, see 



Schleicher et all (|2008ah : ISeager et all (|l999t ). 

The parametrized case B recombination coefficient 
for atomic hydrogen an is given by 



an = F X 10 



-13 



at" 



1 



cm'^ s 



(8) 



with a 4.309, b = -0.6166, c = 0.6703, d = 
0.5300 and t = T/lOf K, which is a fit given by 



Peauignot et al.l (|l99lh . This coefficient takes into 
account that direct recombination into the ground 
state does not lead to a net increase in the num- 
ber of neutral hydrogen atoms, since the photon 
emitted in the recombination process can ionize 
other hydrogen atoms in the neighbourhood. The 
fudge factor F = 1.14 serves to speed up recom- 
bination and is determined from comparison with 
the multilevel code. 

The chemical evolution of the primordial gas 
is solved with a system of rate equations for H~, 



H+, H2, HeH+, D, D+, D", HD+and HD, which 
is largely based on the reaction rates presented 
bv lSchleicher et aL 1 lioosb). For the mutual neu- 
tralization rate of H~ and H+ , we u se the more 
recent result of Stenrup et al. ( 20091 ) . We note 



that the collisional dissociation rates of H2 in this 
compilation are so mewhat larger than those of 
Standi et al. ( 1998h . For this reason, the final H2 
abundance in our calculation is smaller by up to 
an order of magn itude than in the calculation of 
Sethi et al. I j2om which was based on the rates 



of Standi et al.l ( 19981 ). The evolution of the mag- 



netic field strength can be determined from the 
magnetic field energy Eg ~ B'^ /8tt and is followed 
as 

dEe _4:dpEB 



(9) 



dt 3 at p 

The first term describes the evolution of the mag- 
netic field in a homogeneous universe in the ab- 
sence of specific magnetic energy generation or dis- 
sipation mechanisms. If this term dominates, then 
the magnetic field strength evolves as i? oc p^^^ 
due to the expansion of space, and the magnetic 
energy scales with redshift as (1 -|- z)"^, the same 
scaling as for the radiation energy density. The 
second term accounts for corrections due to en- 
ergy dissipation via AD. Such energy dissipation 
is reflected in the evolution of the magnetic Jeans 
mass, which in general is defined as 



Mf 



27r%/G' 



(10) 



with the magnetic flux $ = nr^B, where G is New- 
ton's constant, and an appropriate length scale r. 
For the IGM, a characteristic Jeans length can be 
derived from the Alfven velocity and the dynam- 
ical timescale, yielding ( Subrarnanian &: Barrow 



19981 : ISethi fc Subramanianll2005l) 



Mj 



10^"Mo 



Bn 



3 nG 



(11) 



This equation is formally independent of density, 
as it is expressed here in terms of comoving quan- 
tities. The constant comoving IGM density is thus 
absorbed in the overall normalization. The addi- 
tional heat input affects the evolution of the ther- 
mal Jeans mass as well, which we evaluate as 



Mj 



{7) 



-1/2 



f 5kBT 
\2ijGmp 



3/2 



(12) 
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Here, ks denotes Boltzmann's constant, fi the 
mean molecular weight and rup the proton mass. 

To describe virialization in the first minihalos, 
we employ the spherical collapse model for pres- 
sureless dark matter until an overdensity of 200 
is reached. For this purpose, we first calculate cos- 
mic time in the universe as 



tr. 



:l0g 



1 



cos ( 



sin 9 



(13) 



where is the Hubble constant at z = 0, J^a 
the cosmological density parameter for vacuum en- 
ergy, and is defined via 



tan 6 ~ 



(1 + ^) 



3/2 



(14) 



On the other hand, in the spherical collaps e model, 
we ha ve the following system of equations ( Peebles! 
19931) : 



R 

t 

Al 



Ai{l — cos?]), 
^2(77 - sin??), 
GMAl. 



(15) 
(16) 
(17) 



Here, R is the radius of the collapsing cloud, t 
is the time parameter, Ai and A2 are parameters 
that normalize the length and time scales in this 
problem, M the total mass, which we choose as 
M = 10'^ Mq. (The following results hardly de- 
pend on this choice). With A2 = t2o/27r, where 
^20 is the age of the universe at 2: = 20, we en- 
sure that the cloud collapses at z = 20. Eq. ((T7)) 
then fixes the remaining parameter Ai. Equating 
Eq. p3p and Eq. (fT6|l allows one to derive the pa- 
rameter rj and to calculate the overdensity p/ Pb in 
the protocloud. The latter is given as 



P_ 

Pb 



9 (77 — sin 77)^ 
2(1 — cos 77)'^ 



sin 77) 



^ 9 sin 7/(7; 

2 (1 — cos7y)2 



(18) 

We assume that the formation of the protocloud 
may affect the coherence length Ls of the mag- 
netic field if it is frozen into the gas. We therefore 
adopt Lb = min(i?, k:^^^{\ + z)~^) at this evolu- 
tionary stage. 

2.2. Results 

With the model described above, we calculate 
the evolution of the temperature and chemistry of 
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Fig. 1. — The evolution of the comoving magnetic 
field strength due to AD as a function of redshift 
for different initial comoving field strengths, from 
the homogeneous medium at z = 1300 to virial- 
ization at z = 20. 

the IGM from z = 1300 until virialization in the 
first halos at z = 20 for different initial comov- 
ing field strengths. We note that in our model, 
the comoving field strength is a function of time, 
as the magnetic field is dissipated through AD. In 
the figures given here, Bq therefore labels the ini- 
tial comoving field strength used to initialize the 
calculation. 

As shown in Fig. ([1]) , ambipolar diffusion is par- 
ticularly important for fields with initial comoving 
field strengths of 0.2 nG or less. For stronger fields, 
the dissipation of only a small fraction of their en- 
ergy increases the temperature and the ionization 
fraction of the IGM to such an extent that AD 
becomes less effective. The evolution of temper- 
ature and ionization degree is given in Figs. ([2]) 
and ([3]). While for comoving field strengths up 
to ^ 0.1 nG, the additional heat from ambipolar 
diffusion is rather modest and the gas in the IGM 
still cools below the CMB temperature due to adi- 
abatic expansion, it can increase significantly for 
stronger fields and reaches ~ 10"* K for a comov- 
ing field of 1 nG. At that temperature, Lyman 
a cooling is efficient, and coUisional ionization in- 
creases the ionization degree, which makes AD less 
effective and prevents a further increase in tem- 
perature. Fig. ([3]) also shows the evolution of the 
H2 and HD abundances. As a result of the en- 



6 




10 I , , , ^ 

1000 100 
Redshift 



Fig. 2. — The gas temperature evolution in the 
IGM as a function of redshift for different comov- 
ing field strengths, from the homogeneous medium 
at z = 1300 to virialization at z = 20. For the case 
with Bo = 0.01 nG, we find no difference in the 
thermal evolution compared to the zero-field case. 

hanced ionization fraction, they are increased for 
non-zero field strengths at the point of virializa- 
tion. For relatively weak fields, the abundances 
are increased at all redshifts z < 700 with respect 
to the zero-field case. For field strengths above 
~ 0.7 nG, the heating is so strong that coUisional 
dissociation of H~ and H2 is very efficient and the 
abundances of H2 and HD may even drop below 
the zero-field value. After this point, however, AD 
is shut down by the increased ionization fraction 
and the gas temperature decreases somewhat. The 
molecules thus reform and the abundances again 
become larger than in the zero-field case. 

The evolution of the comoving field strength 
is also reflected in the magnetic Jeans mass, as 
can be seen in Fig. ([4]). Depending on the field 
strength, Mj may become as large as ~ 10^ Mq, 
thus suppressing star formation in the first mini- 
halos. Due to the steep scaling with the comoving 
field to the third power, it decreases quickly if the 
magnetic field is weaker. Also the thermal Jeans 
mass is affected in the presence of magnetic fields 
due to the increase of the gas temperature. The 
impact on Mj is plotted in Fig. ([5]). The increase 
is significant for comoving field strengths above 
0.1 nG. Up to ~ 1 nG, the thermal Jeans mass 
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Fig. 3. — The evolution of ionization degree, H2 
and HD abundances as a function of redshift for 
different comoving field strengths, from the homo- 
geneous medium at z = 1300 to virialization at 
z — 20. For the case with Bq = 0.01 nG, we find 
no difference in the chemical evolution compared 
to the zero-field case. 

dominates over the magnetic one. For stronger 
fields, the gas temperature cannot increase much 
further, as AD is shut down, while the magnetic 
Jeans mass would still scale with Bq and would 
therefore dominate. 

In summary, our results show that the presence 
of magnetic fields may affect the chemical initial 
conditions for star formation in the first halos at 
z ~ 20. They further indicate that the magnetic 
and/or the thermal Jeans mass may be consid- 
erably increased for comoving fields of at least 
0.1 nG, implying that small halos need to accrete 
more mass until they can eventually collapse. 

3. The protostellar collapse phase 

In this section, we use the results from the pre- 
vious section, in particular the chemical initial 
conditions and the dissipation of the magnetic en- 
ergy, as input parameters to calculate the chem- 
istry during the collapse of the first protostars. 
This calculation therefore assumes that the corre- 
sponding halos are massive enough, such that the 
gas is gravitationally unstable and can form stars. 
If the halo is less massive initially, the gas would 
linger at low densities until enough material is ac- 



7 




Fig. 4. — The evolution of the magnetic Jeans 
mass as a function of redshift for different comov- 
ing field strengths, from the homogeneous medium 
ai z = 1300 to virialization at z = 30. 

creted to go into collapse. Such accretion could be 
due to minor or major mergers as well as cooling 
flow activity. 

3.1. Model assumptions 

To follow the evolution of the primordial gas 
during the protostellar collapse phase, we em- 
plo y and extend t h e one -zone model developed 



plo y ana extena t ne one-zone model aevelopea 
bv iGlover &: SavinI Hooi). This model assumes 



that the collapse occurs on the free-fall timescale, 
unaffected by changes in the thermal energy of 
the gas. It incorporates the most extensive treat- 
ment of primordial gas chemistry currently avail- 
able, involving 392 reactions amongst 30 atomic 
and molecular species, and hence allows us to 
follow the fractional ionization of the gas accu- 
rately even for values as small as Xc ~ 10^^^. 
The chemical rate equations are evolved simulta- 
neously with the thermal energy equation, and the 
coupled set of equations are solved implicitly using 



couplea set ol equations are solvea implicitly using 
the DVODE integrator terown et all 1 1989( 1. The 
effects of chemical heating due to three-body H2 
formation and chemical cooling due to H2 and HD 
dissociation are included, and play an important 
role in regulating the temperature of the gas at 
densities n > 10* cm~^. Particularly important 
for this application is the fact that it correctly 
models the evolution of the ionization degree and 



Fig. 5. — The evolution of the thermal Jeans mass 
as a function of redshift for different comoving field 
strengths, from the homogeneous medium at z = 
1300 to virialization at z = 20. For the case with 
Bq = 0.01 nG, we find no difference in the thermal 
Jeans mass compared to the zero-field case. 

the transition at densities of ~ 10* cm~'^ where 
Li"*" becomes the main charge car rier. Further de- 
tails o f the model can be found in lGlover fc Savin 
(I2OO9I) . 

The additional heat input due to magnetic en- 
ergy dissipation may however have an important 
impact on the thermal evolution of the gas and 
delay gravitational collapse. It may thus lead to a 
backreaction on the heating rate for gravitational 
contraction, as for a given change in volume, the 
energy input from pdV work is fixed, while the 
amount of energy that is radiated away increases if 
collapse is delayed. At the same time, the amount 
of energy dissipated by AD may be increased. To 
model the complex interplay of the heating and 
cooling rates involved, we therefore need to eval- 
uate how the collapse timescale is affected by the 
thermodynamics of th e gas. For this purpo se, we 
adopt the approach of lOmukai et al. based 
on the Larson-Penston type self-similar solution 
(iLarson Il969l: iPenstonI 1 19691 ) as generalized by 
lYahill (|1983^ In this solution, the actual collapse 
timescale icoii is related to the free-fall timescale 
tg as 



icoU — 



1 



(19) 
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where / describes the ratio between the pressure 
gradient force and the gravitational force at the 
center. It can be calculated from an effective equa- 
tion of state parameter 7 = d\ogp/d\og p as 



/ 



0, 7 < 0.83, 



(20) 



= 0.6 + 2.5(7-1)- 6.0(7 0.83 <7<1, 
= 1.0 + 0.2(7 -4/3)- 2.9(7 -4/3)^ 7>1- 



At 7 = 4/3, / reaches unity and the collapse 
timescale diverges, indicating that the self-similar 
solution is no longer valid under these circum- 
stances. In reality, a hydrostatic core would form 
and thereaf ter contract during i ts further evolu- 
tion. As in lOmukai et all (|2005f) . we mimic this 
effect by adopting an upper limit / = 0.95. 

During protostellar collapse, magnetic fields are 
typically found to scale as a power-law with den- 
sity p. Assuming ideal MHD with flux freezing 
and spherical collapse, one expects a scaling with 
/O^/^ in the case of weak fields. Deviations from 
spherical symmetry such as expected for dynam- 
ically important fields give rise to shallower sca l- 



ings, e.g. B oc 0°; ^ (Baneriee fc Pudrit d I2OO6I) 



B PC o^/^ (Spitzer 1981 : iHcnncbe lle fc Fromanel 
20081: (Hennebelle fc Tevssier. ,2008) . This is be- 
cause collapse preferentially occurs along the field 
lines in the latter case, and is slowed down in the 
perpendicular direction. 

The difference in these scaling behaviors may 
seem negligible at first sight, and is indeed hard 
to distinguish in simulations that often cover only 
a few orders of magnitude in density. For our 
calculations, however, it is an important issue, as 
they cover about 12 orders of magnitude i n den- 



sity. As the simulations of Machida et al. ( 20061 ) 



cover a similar range of densities and magnetic 
field strengths as our calculations, we use them to 
derive an appropriate scaling relation. Of course, 
this relation is only approximate, as their simu- 
lation does not include a consistent treatment of 
non-equilibrium chemistry and the additional heat 
input from AD, which may have a backreaction on 
the dynamics. In their results, the slope a of the 
scaling relation B (x p" depends on the ratio of the 
thermal Jeans mass to the magnetic Jeans mass, 
which we calculate from Eqs. (fTO|) and (fT2)) . We 
find the empirical scaling law 



0.57 



Mj 
Mf 



0.0116 



(21) 



This relation matches our expectations described 
above. For negligible magnetic pressure and 
low magnetic Jeans masses, the magnetic field 
strength increases more rapidly with density, while 
for strong magnetic pressure, gas collapse occurs 
preferentially along the magnetic field lines, and 
the scaling relation flattens. Even though the rela- 
tion depends only weakly on this ratio, one should 
note that Mj/Mf may vary by several orders of 
magnitude depending on the initial field strength. 
As an additional caveat, we point out that this 
relation may to some extent depend on the initial 
conditions and change depending on the individual 
properties of the first minihalos. For this reason, 
we will not only explore the consequences of this 
particular relation, but also consider the effects 
from softer or steeper relations below. 

The calculation of the magnetic Jeans mass as 
given in Eq. (|10p requires an assumption regard- 
ing the length scale r of the dense region. Nu- 
merical hydrodynamics simulations show that re- 
gions of a given density have length scales com- 



parablc to the th ermal Jeans length (|Abel et al 



i2002 : Bromm fc La rson 2004) . We expect similar 
effects from the magnetic Jeans mass and there- 
fore take the maximum of the thermal and mag- 
netic Jeans length. As shown in the previous sec- 
tion, the thermal Jeans mass dominates over the 
magnetic Jeans mass for comoving field strengths 
up to '--^ 1 nG. We therefore initialize this length 
scale with the thermal Jeans length at the begin- 
ning of collapse, and follow its further evolution 
by taking the maximum of the thermal and the 
magnetic Jeans length. 

At every timestep, the magnetic field strength 
is updated as 



Br. 



B, 



Id 



Pn 



Pold 



(22) 



The coefhcient a is calculated at every timestep 
such that a change in the ratio of the thermal 
and magnetic Jeans mass, either due to chemistry, 
magnetic energy dissipation or magnetic field am- 
plification, may lead to a backreaction on the in- 
crease of the field strength with density. However, 
Eq. ([221) does not incorporate the effect of mag- 
netic energy dissipation on the field strength. As 
for the IGM, we can calculate the AD heating rate 
from Eq. As we are here considering a much 
larger range of densities, additional processes need 
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to be taken into account to calculate the AD re- 
sistivity correctly. In particular, at a density of 
~ 10® cni~^, the three-body H2 formation rates 
start to increase the H2 abundance significantly, 
such that the gas is fully molecular at densities 
of ~ 10^^ cm~^. As a further complication, the 
proton abundance drops considerably at densities 
of ~ 10® cm~^, such that Li+ becomes the main 



10000 



charge carrier (jGlover fc Savinll2009l ). To model 
the total AD resistivity correctly, we therefore 
need to calculate the resistivities of H, He and H2 
and take into account collisions with protons and 
Li+. For collisions with proton s, we adopt the mo - 
mentum transfer coefficients of lPinto et al.l (j2008l ). 
while we use the polar ization approximati on for 
collisions with Li+ fsee lPinto k Gallill2008h . 

With these considerations, we can calculate the 
AD heating rate and its impact on the thermody- 
namics. We evaluate its impact on the magnetic 
field strength as 



Bear = VStT^ Bl^^^J^-K - SILaD, 



(23) 



where St denotes the timestep, and -Bcor and 
-Buncor denote the magnetic field corrected and 
not corrected for the effects of AD, respectively. 
The equation follows from the dependence of the 
magnetic energy density on the field strength and 
the impact of AD on the magnetic energy den- 
sity. Finally, we checked whether Ohmic diffusion 
(see lPinto et al.l 120081) can be relevant under pri- 
mordial conditions, finding that AD is always the 
dominant source of energy dissipation. 

3.2. Results 

With the model described above, we calculate 
the evolution of the magnetic field strength and 
its impact on protostellar collapse in primordial 
gas. For this purpose, we adopt the results from 
§[2] for the field strength, the temperature and the 
chemical abundances at virialization and use them 
as initial conditions for the collapse calculation. 
For comoving field strengths of at least 0.1 nG, 
we showed in the previous section that the abun- 
dances of free electrons, H2 and HD may be sig- 
nificantly increased in the IGM. During collapse, 
the enhanced fraction of free electrons may fur- 
ther catalyze the formation of molecules, so that 
the cooling rate increases in the presence of mag- 
netic fields. On the other hand, there is also more 
heat input from magnetic energy dissipation. In 
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Fig. 6. — The gas temperature as a function of 
density for different comoving field strengths. For 
Bq = 0.01 nG, the thermal evolution corresponds 
to the zero-field case. 



the absence of magnetic fields, it was shown that 
such initial condi tions may lead to coo ling down to 
the CMB floor (lYoshida et aLlbOOTaf bl) and char- 
acteristic stellar masses of ~ 10 Mq. We will as- 
sess here how these results change in the presence 
of magnetic energy dissipation. Finally, note that 
the label Bq refers to the comoving field strength 
used to initialize the IGM calculation. The corre- 
sponding physical field strength at the beginning 
of the collapse phase is given in Table [T] 



Bq [nG] B [nG] 

i 4.8 X 10^ 

0.3 6.3 X 102 

0.1 1.3 xlO^ 

0.01 1.0 X 10^ 



Table 1: The physical field strength B at begin- 
ning of collapse as a function of the comoving field 
strength Bq used to initialize the IGM calculation 
at z = 1300. 

Fig. [5] shows the temperature evolution as a 
function of density for different comoving field 
strengths. For comoving fields of 0.01 nG or less, 
there is virtually no difference in the temperature 
evolution from the zero- field case. For comoving 
fields of ~ 0.1 nG, cooling wins over the additional 
heat input in the early phase of collapse, and the 
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Fig. 7. — Ionization degree, H2 and HD abundance 
as a function of density for different comoving field 
strengths. For Bq — 0.01 nG, the thermal evolu- 
tion corresponds to the zero- field case. 



temperature decreases slightly below the zero-field 
value at densities of 10'^ cm~'^. At higher densities, 
the additional heat input dominates over cooling 
and the temperature steadily increases. At den- 
sities of ~ 10^ cm~'^, the abundance of protons 
drops considerably and increases the AD resistiv- 
ity defined in Eq. ([3]) and the heating rate until Li+ 
becomes the main charge carrier. In particular for 
comoving fields larger than ~ 0.1 nG, this transi- 
tion is reflected by a small bump in the tempera- 
ture evolution due to the increased heating rate in 
this density range. The transition is also visible in 
Fig. [51 which shows the evolution of magnetic field 
strength with density. At the transition, we find a 
flattening of the relation between density and field 
strength, as a significant fraction of the energy can 
be dissipated at this stage. As expected from the 
evolution of the H2 abundance, the main contribu- 
tion to the total resistivity is due to the resistivity 
of atomic hydrogen, until the gas becomes fully 
molecular at densities of ~ 10"'^^ cm~^. 

The evolution of the abundances of free elec- 
trons, II2 and HD are given in Fig.[71 With higher 
initial field strength, the abundance of free elec- 
trons increases, inducing the formation of more 
molecules. The molecular abundances thus in- 
crease with field strength as well. Interestingly, 
the molecular fraction becomes independent of 



Fig. 8. — The physical field strength as a function 
of density for different comoving field strengths. 



field strength at densities of ^ 10 cm~ where 
three-body II2 formation takes over. Note that 
there are still small differences in the electron 
abundance. 

Apart from the transition where Li"*" becomes 
the dominant charge carrier, the magnetic field 
strength usually increases more rapidly than p"'^, 
and weak fields increase more rapidly than strong 
fields. This is what one naively expects from 
Eq. (|2ip . and it is not significantly affected by 
magnetic energy dissipation. Another important 
point is that comoving fields of only 10^^ nG 
are amplified to values of ~ 1 nG at a density 
of 10^ cm~^. Such fields are required to drive 
protostellar outfiows that can magnetize the IGM 
IjMachida eral1l2006h . 

Fig. [5] shows the evolution of the thermal and 
magnetic Jeans mass during collapse. The ther- 
mal Jeans masses are quite different initially, but 
as the temperatures reach the same order of mag- 
nitude during collapse, the same holds for the ther- 
mal Jeans mass. The thermal Jeans mass in this 
late phase has only a weak dependence on the 
field strength. As expected, the magnetic Jeans 
masses are much more sensitive to the comoving 
field strength, and initially differ by about two or- 
ders of magnitude for one order of magnitude dif- 
ference in the field strength. For comoving fields 
of ~ 1 nG, the magnetic Jeans mass dominates 
over the thermal one and thus determines the mass 
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Fig. 9. — Thermal (thin hnes) and magnetic (thick 
Hnes) Jeans mass as a function of density for differ- 
ent comoving field strengths. For Bq = 0.01 nG, 
the thermal evolution and thus the thermal Jeans 
mass corresponds to the zero-field case. 

scale of the protocloud. For ~ 0.3 nG, both masses 
are roughly comparable, while for weaker fields 
the thermal Jeans mass dominates. The magnetic 
Jeans mass shows features both due to magnetic 
energy dissipation, but also due to a change in the 
thermal Jeans mass, which sets the typical length 
scale and thus the magnetic flux in the case that 
Mj > Mf. 

We have checked the sensitivity of our results 
to our model assumptions, finding that the cor- 
rect scaling relation between the magnetic field 
strength and density is crucial for these results. 
As an example, we show how the results for a 1 nG 
field depend on this scaling relation in Figs. PH)) 
and pTjl . In this particular case, the fitting for- 
mula in Eq. (|2ip matches best with a constant 
scaling a = 0.55. For larger values of a, the 
heating rate increases rapidly with density, and 
reaches a critical value for which the H2 dissocia- 
tion rates become very effective. These dissocia- 
tion rates change by several orders of magnitude in 
this temperature range are extremely sensitive to 
the temperature evolution. Once activated, they 
act to reduce the cooling from H2, which amplifies 
the temperature increase further and leads to the 
sudden jump shown in Fig. (|10p . On the other 
hand, for a < 0.55, the temperature might be 
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Fig. 10. — The gas temperature as a function of 
density for a comoving field strength of 1 nG and 
different assumptions for the scaling between den- 
sity and field strength, both assuming constant 
slopes a and our fiducial model based on the fit 
given in Eq. ((2T|) . 

smaller than in our standard runs. It is interest- 
ing to note that the bump due to the transition 
where Li"*" becomes the main charge carrier be- 
comes less prominent for steeper relations, as the 
latter considerably enhances the heating rate at 
higher densities. 

We have also checked the sensitivity to the scal- 
ing relation for a comoving field strength of 0.1 nG 
(see Fig. [T^ . In this case, the differences are 
smaller, but still significant at densities larger than 
10® cm^"^. For this initial field strength, the bump 
at 10^ cm"'^ becomes more prominent for a steeper 
relation. This is because AD heating is only sig- 
nificant for a > 0.55. For smaller values of a, the 
differences with the zero-field case are mostly due 
to the change in the initial chemical abundances, 
and not due to the effects of AD during the col- 
lapse. 

4. Discussion 

The calculations above show that comoving 
magnetic fields of ^ 0.1 nG or more significantly 
increase the thermal and magnetic Jeans mass in 
the IGM to 10^ - 10^ Mq. The thermal evolution 
is modified both due to changes in the chemical 
initial conditions at the beginning of the collapse 
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Fig. 11. — Ionization degree (bold lines) and H2 
abundance (thin lines) as a function of density for 
a comoving field strength of 1 nG and different 
assumptions for the scaling between density and 
field strength, both assuming constant slopes a 
and our fiducial model based on the fit given in 
Eq. (Ell). 



phase, as well as due to the further energy input 
from magnetic energy dissipation. A correct mod- 
eling of the AD resistivities is quite important, 
as these rise considerably at high densities when 
the proton abundance drops before Li+ becomes 
the main charge carrier. Our results further in- 
dicate that comoving field strengths of 10~^ nG 
are sufficiently amplified during protostellar col- 
la pse to reach th e criti cal field strength derived 
bv lMachida et al. for the formation of jets, 

which may magnetize the IGM. Of course, addi- 
tional and perhaps even more important contri- 
butions may be generated once massive Pop. Ill 
stars are formed ( Bisnovatvi-Kogan et al.il973al bi: 



Bisnovatvi-Kogan fc Lamzin 19771 ) 

During protostellar collapse, we do not find a 
significant temperature change at the minimum 
near densities of ~ 10'^ cm~'^, and therefore do not 
expect a large change in the characteristic frag- 
mentation mass scale. At higher densities, the ad- 
ditional heat input will slow down the collapse, in 
particular at densities of ~ 10^ cm~'^ where Li+ 
becomes the main charge carrier. This may lead to 
deviations from spherical symmetry and favor the 
formation of flattened structures like an accretion 



Fig. 12. — The gas temperature as a function of 
density for a comoving field strength of 0.1 nG 
and different assumptions for the scaling between 
density and field strength, both assuming constant 
slopes a and our fiducial model based on the fit 
given in Eq. ((2T|) . 



disk. Due to the higher temperatures, a higher ac- 
cretion rate onto the protostar may be expected. 
Part of this may however be balanced by outflows 
induced by the magnetic fields, and the dynamical 
implications of our results need to be explored in 
three-dimensional MHD simulations. 

In summary, our results confirm that the pres- 
ence of magnetic fields may increase the Jeans 
mass in the IGM considerably, which may sup- 
press gravitational collapse in the first minihalos 
and therefore delay re ionization. As shown by 
Schleicher et al.l (|2008al ). it is therefore possible to 
derive upper limits on the comoving field strength, 
which are of the order of a few nG, depending on 
assumptions for the stellar population. With sub- 
sequent works, we plan to further reduce the un- 
certainties regarding the stellar population, such 
that tighter constraints can be derived. In addi- 
tion, we expect that the Planck data and its more 
accurate measurement of the reionization optical 
depth will be very valuable in this respect. 

It is also of interest to understand how the ef- 
fects of magnetic fields depend on the environmen- 
tal conditions. For example, in relic HII regions 
and atomic cooling halos, the temperature and 
ionization fraction could be significantly increased, 
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leading to a difFerent evolution of the chemistry 



(jYoshida et al.M2007al |hl). We have checked that 
for comoving field strength of ~ 1 nG, the addi- 
tional heat input from magnetic energy dissipation 
dominates over the additional cooling from the en- 
hanced H2 and HD abundance, while for comoving 
fields of < 0.1 nG, the difference from the zero-field 
case is small, as the higher ionization degree makes 
AD less effective. Further effects may be due to 
the presence o f turbulence in these massive atomi c 
coohng halos (jClark et al.ll2008HGreif et al.ll2008f) . 
In particular, the coherence length of the magnetic 
field may be reduced increasing the AD heating 
rate, while energy dissipation from decaying MHD 
turbulence may deposit further heat into the gas. 

Subsequent generations of stars may form 
from gas that has been enriched with metals, 
and presumably stronger magnetic fields due to 
stellar winds or supernova explosions. The im- 
pact of the additional metallicity ha s been dis 



cussed extensively in the literature (.Bromm et al 
20011 ISchneider et all |2005 lOmukai et all 1 2005 



Schneider et al." 2006'; Glover & JapDseri '20071 
JaPDScn ct al. 2007, 2009a, b: .Tornatore et al. | 
2007 [IClark et al.l2008HOmukai et al ll2008HSmith et al 
20091) . 



Conversely, both the amount of magnetic 
energy generated in the IGM by the first gen- 
eration of stars, as well as their possible effects 
on the subsequent generatio ns a r e still largely 



unexplored (but see Rees 1987t Kandus et al 
I2OO4 ). At low metallicity, d ust cooling becomes 



important at high densities (jOmukai et al. I I2OO5I: 



[Schncide r ct al. .2006. ) and may easily dominate 
over the additional heating from ambipolar dif- 
fusion. At lower densities, the main effect of 
dust is to increase the abundance of molecular 
hydrogen due to H2 fo rmation on dust grains 
( Cazaux fc Spaaiiil2009l ). This may also increase 
the cooling rate significantly. At the same time, 
the presence of metallic ions may increase the ion- 
ization fraction in the gas, such that ambipolar 
diffusion may become less important. However, 
it also implies that the dynamical effects of mag- 
netic fields will play a more important role, as the 
coupling between the gas and the magnetic field 
will be more efficient. 

Finally, we stress a fundamental difference in 
the magnetically controlled collapse of clouds of 
primordial and non-primordial chemical compo- 
sitions. Whereas the electric charge at densities 



n > 10^ cm""^ in a primordial cloud is mostly car- 
ried by Li+ ions and electrons (see Sect. 13.21 ). in a 
cloud of solar chemical composition the charged 
component in the same density range is domi- 
nated by negatively and positively charge d dust 



grain s (jDesch fc Mouschovi as 2001; Nakan o et al 
20021 ). This has important consequences on the 



evolution of the magnetic field because the elec- 
tric resistivities associated to ambipolar diffusion. 
Hall effect, and Ohmic dissipation are expected 
to differ even by orders of magnitud e in the two 
cases, as shown e.g. by IPinto et"al ] (I2OO8) . Un- 
fortunately, while the polarization approximation 
adopted in our calculation of the AD resistivity ap- 
pears to be justified for collisions of Li"*" with He 



( Cassidv k ElfordI 119851) a nd H2 (jRoeggen et al 



2OO2I : [Dickinson et al.lll982D . collisions of Li+ with 
H atoms, especially relevant for the primordial gas, 
have not been studied in detail cither theoretically 
or experimentally. 

As noted above, all of these results are sensi- 
tive to the correct scaling relation between the 
magnetic field strength and the gas density. A 
correct and self-consistent derivation of this rela- 
tion is therefore important for our understanding 
of primordial star formation. 

D.R.G.S. thanks the Heidelberg Graduate 
School of Fundamental Physics (HGSFP) and the 
State of Baden- Wiirttemberg for financial sup- 
port. R.S.K. acknowledges support from the Ger- 
man Science Foundation (DFG) under the Emmy 
Noether grant KL1358/1 and the Priority Pro- 
gram SFB 439 Galaxies in the Early Universe. 
S.C.O.G. acknowledges funding from the DFG via 
grant KL1358/4. RB is funded by the DFG under 
the grant BA 3706/1. This work was supported in 
part by a FRONTIER grant of Heidelberg Univer- 
sity sponsored by the German Excellence Initia- 
tive. We also acknowledge funding from the Eu- 
ropean Commission Sixth Framework Programme 
Marie Curie Research Training Network CON- 
STELLATION (MRTN-CT-2006-035890). We 
thank the anonymous referee for many clarifying 
remarks that helped to improve this paper. 

REFERENCES 

Abel, T., Anninos, P., Zhang, Y., & Norman, 
M. L. 1997, New Astronomy, 2, 181 



14 



Abel, T., Bryan, G. L., & Norman, M. L. 2002, 
Science, 295, 93 

Anninos, P., Zhang, Y., Abel, T., & Norman, 
M. L. 1997, New Astronomy, 2, 209 

Bamba, K., Ohta, N., & Tsujikawa, S. 2008, 
Phys. Rev. D, 78, 043524 

Banerjee, R., & Jedamzik, K. 2004, Phys. Rev. D, 
70, 123003 

Banerjee, R., & Pudritz, R. E. 2006, ApJ, 641, 949 

Banerjee, R., & Pudritz, R. E. 2007, ApJ, 660, 479 

Baym, G., Bodeker, D., & McLerran, L. 1996, 
Phys. Rev. D, 53, 662 

Beck, R. 2009, Ap&SS, 320, 77 

Berkhuijsen, E. M., Beck, R., & Hoernes, P. 2003, 
A&A, 398, 937 

Bernet, M. L., Miniati, F., Lilly, S. J., Kronbcrg, 
P. P., & Dessauges-Zavadsky, M. 2008, Nature, 
454, 302 

Bertolami, O., & Mota, D. F. 1999, Physics Let- 
ters B, 455, 96 

Biermann, L. 1950, Zeitschrift Naturforschung 
Teil A, 5, 65 

Bisnovatyi-Kogan, G. S., Ruzmaikin, A. A., & 
Syunyaev, R. A. 1973a, Soviet Astronomy, 17, 
137 

— . 1973b, AZh, 50, 210 

Bisnovatyj-Kogan, G. S., & Lamzin, S. A. 1977, in 
Early Stages of Stellar Evolution, ed. K. Kaerre 
& H. Lindkvist, 107-118 

Bouwens, R. J., & lUingworth, G. D. 2006, Nature, 

443, 189 

Bromm, V., Ferrara, A., Coppi, P. S., & Larson, 
R. B. 2001, MNRAS, 328, 969 

Bromm, V., & Larson, R. B. 2004, ARA&A, 42, 
79 

Brown, P. N., Byrne, G. D., & Hindmarsh, A. C. 
1989, SIAM J. Sci. Stat. Comput., 10, 1038 

CampaneUi, L. 2008, ArXiv 0805.0575 



Campanelli, L., Cea, P., Fogli, G. L., & Tedesco, 
L. 2008, Phys. Rev. D, 77, 043001 

Cassidy, R. A., & Elford, M. T. 1985, Australian 
Journal of Physics, 38, 587 

Cazaux, S., & Spaans, M. 2009, A&A, 496, 365 

Cheng, B., & Olinto, A. V. 1994, Phys. Rev. D, 
50, 2421 

Ciardi, B., & Ferrara, A. 2005, Space Science Re- 
views, 116, 625 

Clark, P. C, Glover, S. C. O., & Klessen, R. S. 
2008, ApJ, 672, 757 

Desch, S. J., & Mouschovias, T. C. 2001, ApJ, 550, 
314 

Dickinson, A. S., Lee, M. S., & Lester, Jr., W. A. 
1982, Journal of Physics B Atomic Molecular 
Physics, 15, 1371 

Dolag, K., Bartelmann, M., & Lesch, H. 2002, 
A&A, 387, 383 

Dolag, K., Grasso, D., Springel, V., & Tkachcv, I. 
2005, Journal of Cosmology and Astro-Particle 
Physics, 1, 9 

Glover, S. 2005, Space Science Reviews, 117, 445 

Glover, S. C. O., & Jappsen, A.-K. 2007, ApJ, 666, 
1 

Glover, S. C. O., & Savin, D. W. 2009, MNRAS, 
393, 911 

Grasso, D., & Rubinstein, H. R. 1996, Physics Let- 
ters B, 379, 73 

— . 2001, Phys. Rep., 348, 163 

Grcif, T. H., Johnson, J. L., Klessen, R. S., & 
Bromm, V. 2008, MNRAS, 387, 1021 

Han, J. L. 2008, Nuclear Physics B Proceedings 
Supplements, 175, 62, Proceedings of the XIV 
International Symposium on Very High Energy 
Cosmic Ray Interactions, Weihai, China 

Hanasz, M., Otmianowska-Mazur, K., Kowal, G., 
& Lesch, H. 2009, A&A, 498, 335 

Hennebelle, P., & Fromang, S. 2008, A&A, 477, 9 



15 



Hennebelle, P., & Teyssier, R. 2008, A&A, 477, 25 

lye, M., Ota, K., Kashikawa, N., Furusawa, H., 
Hashimoto, T., Hattori, T., Matsuda, Y., Mo- 
rokuma, T., Ouchi, M., & Shimasaku, K. 2006, 
Nature, 443, 186 

Jappsen, A.-K., Glover, S. C. O., Klessen, R. S., 
& Mac Low, M.-M. 2007, ApJ, 660, 1332 

Jappsen, A.-K., Klessen, R. S., Glover, S. C. O., 
& Mac Low, M.-M. 2009a, ApJ, 696, 1065 

Jappsen, A.-K., Low, M.-M. M., Glover, S. C. O., 
Klessen, R. S., & Kitsionas, S. 2009b, ApJ, 694, 
1161 

Jedamzik, K., Katalinic, V., & Olinto, A. V. 1998, 
Phys. Rev. D, 57, 3264 

Kandus, A., Opher, R., & Barros, S. R. B. 2004, 

Braz. J. Phys., 34, 4b 

Kim, E.-J., Olinto, A. V., & Rosner, R. 1996, ApJ, 
468, 28 

Kronberg, P. P. 1994, Reports on Progress in 
Physics, 57, 325 

Kulsrud, R., Cowley, S. C., Gruzinov, A. V., & 
Sudan, R. N. 1997, Phys. Rep., 283, 213 

Kulsrud, R. M., & Zweibel, E. G. 2008, Reports 
on Progress in Physics, 71, 046901 

Larson, R. B. 1969, MNRAS, 145, 271 

Machida, M. N., Omukai, K., Matsumoto, T., & 
Imitsuka, S.-I. 2006, ApJ, 647, LI 

Maki, H., & Susa, H. 2004, ApJ, 609, 467 

— . 2007, PASJ, 59, 787 

Nakano, T., Nishi, R., & Umebayashi, T. 2002, 

ApJ, 573, 199 

Omukai, K., Schneider, R., & Haiman, Z. 2008, 
ApJ, 686, 801 

Omukai, K., Tsuribe, T., Schneider, R., fc Ferrara, 
A. 2005, ApJ, 626, 627 

Parker, E. N. 1992, ApJ, 401, 137 

Parker, L. 1968, Physical Review Letters, 21, 562 



Peebles, P. J. E. 1993, Principles of physical cos- 
mology (Princeton Series in Physics, Princeton, 
NJ: Princeton University Press) 

Penston, M. V. 1969, MNRAS, 144, 425 

Pequignot, D., Petitjean, P., & Boisson, C. 1991, 
A&A, 251, 680 

Pinto, C., & Gain, D. 2008, A&A, 484, 17 

Pinto, C., Gain, D., & Bacciotti, F. 2008, A&A, 
484, 1 

Price, D. J., & Bate, M. R. 2008, MNRAS, 385, 
1820 

Quashnock, J. M., Loeb, A., & Spergel, D. N. 

1989, ApJ, 344, L49 

Rees, M. J. 1987, QJRAS, 28, 197 

Riotto, A., & Trodden, M. 1999, Annual Review 
of Nuclear and Particle Science, 49, 35 

R0eggen, I., SkuUerud, H. R., L0vaas, T. H., 
& Dysthe, D. K. 2002, Journal of Physics B 
Atomic Molecular Physics, 35, 1707 

Schleicher, D. R. G., Banerjee, R., & Klessen, R. S. 
2008a, Phys. Rev. D, 78, 083005 

— . 2009, ApJ, 692, 236 

Schleicher, D. R. G., Galli, D., Palla, F., Camen- 
zind, M., Klessen, R. S., Bartelmann, M., & 
Glover, S. C. O. 2008b, A&A, 490, 521 

Schneider, R., Ferrara, A., Salvaterra, R., Omukai, 
K., & Bromm, V. 2003, Nature, 422, 869 

Schneider, R., Salvaterra, R., Ferrara, A., & Cia- 
rdi, B. 2006, MNRAS, 369, 825 

Seager, S., Sasselov, D. D., & Scott, D. 1999, ApJ, 

523, LI 

— . 2000, ApJS, 128, 407 

Seshadri, T. R., & Subramanian, K. 2001, Physi- 
cal Review Letters, 87, 101301 

Sethi, S. K., Nath, B. B., & Subramanian, K. 2008, 
MNRAS, 387, 1589 

Sethi, S. K., & Subramanian, K. 2005, MNRAS, 
356, 778 



16 



Sigl, G., Olinto, A. V., & Jedamzik, K. 1997, 
Phys. Rev. D, 55, 4582 



— . 2008, Science, 321, 669 



Zweibel, E. G. 2006, Astronomische Nachrichten, 
327, 505 



Silk, J., & Langer, M. 2006, MNRAS, 371, 444 



Smith, B. D., Turk, M. J., Sigurdsson, S., O'Shea, 
B. W., & Norman, M. L. 2009, ApJ, 691, 441 

Spitzer, Jr., L. 1981, Physical processes in the in- 
terstellar medium. 

Standi, P. C., Lepp, S., & Dalgarno, A. 1998, ApJ, 
509, 1 

Stenrup, M., Larson, A., & Elander, N. 2009, 
Phys. Rev. A, 79, 012713 

Subramanian, K., & Barrow, J. D. 1998, 
Phys. Rev. D, 58, 083502 

Tan, J. C., & Blackman, E. G. 2004, ApJ, 603, 
401 

Tashiro, H., & Sugiyama, N. 2006, MNRAS, 368, 
965 

Tashiro, H., Sngiyama, N., & Banerjee, R. 2006, 
Phys. Rev. D, 73, 023002 

Tornatore, L., Ferrara, A., & Schneider, R. 2007, 

MNRAS, 382, 945 

Turner, M. S., & Widrow, L. M. 1988, 
Phys. Rev. D, 37, 2743 

Wasserman, I. M. 1978, PhD thesis, AA(Harvard 
University.) 

Wong, W. Y., Moss, A., & Scott, D. 2008, MN- 
RAS, 386, 1023 

Wyithe, J. S. B., & Loeb, A. 2003, ApJ, 588, L69 

Xu, H., O'Shea, B. W., Collins, D. C., Norman, 
M. L., Li, H., & Li, S. 2008, ApJ, 688, L57 

Yahil, A. 1983, ApJ, 265, 1047 

Yamazaki, D. G., Ichiki, K., Kajino, T., & Math- 
ews, G. J. 2006, ApJ, 646, 719 

— . 2008, Phys. Rev. D, 78, 123001 

Yoshida, N., Oh, S. P., Kitayama, T., & Hern- 
quist, L. 2007a, ApJ, 663, 687 

Yoshida, N., Omukai, K., & Hernquist, L. 2007b, This 2-column preprint was prepared with the AAS IATeX 

ApJ, 667, L117 macros v5.2. 



17 



